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Abstract 



The ncutrinoproduction of and Df mesons is studied, using the data obtained with 
the SKAT bubble chamber at the Serpukhov accelerator. It is found that the (p produc- 
tion occurs predominantly in the forward hemisphere of the hadronic c.m.s. (at > 0, 
xp being the Feynman variable), with the mean yield strongly exceeding the expected 
yield of directly produced (j) mesons and varying from {n^{xp > 0)) = (0.92ib0.34)-10^^ 
atW > 2 GeV up to (1.23 ± 0.53) • lO'^ at > 2.6 GeV and (1.44 ± 0.69) • 10"^ at 
W > 2.9 GeV, W being the invariant mass of the hadronic system. The yield of leading 
Df mesons carrying more than z = 0.85 of the current c- quark energy is estimated: 
atW> 2.9 GeV, {n^+{z > 0.85)) = (6.64 ± 1.91) • 10"^. It is shown, that the shape 
of the (p meson differential spectrum on xj? is reproduced by that expected from the 
Df (f)X decays which, however, can account for only the half of the measured (p 
yield. 
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1 Introduction 



The total and differential yields of hadrons in leptonuclear reactions reflect the space-time 
structure of the quark string fragmentation and the formation of hadrons, both produced 
directly or as a result of secondary intranuclear interactions or originated from the decays of 
resonances. The latters play a significant role in the production of stable hadrons (pions and 
kaons). For example, in the case of neutrino-induced reactions the fraction of pions from 
the decays of light meson resonances (up to masses ~ 1 GeV/c^) composes about 20-40%, 
depending on the pion type and the neutrino energy E^, (see [Ij and references therein). The 
information on the space-time structure of the hadron formation based only on the data on 
stable hadrons would, therefore, be rather incomplete and should be complemented by the 
data on the resonance production. The latter data themselves provide more direct infor- 
mation on the quark string fragmentation, since they are, in general, expected to be less 
contaminated by a contribution from the decays of the higher-mass resonances. This contri- 
bution can be rather different for various resonances and in various kinematical regions of 
their production. In particular, in neutrino-induced reactions the production of a favorable 
resonance, which can carry a large fraction of the current quark energy and hence occupies 
the forward hemisphere in the hadronic c.m.s. (with xp > 0, xp being the Feynman vari- 
able), is contributed by the decay of the higher-mass resonances mainly at Xi;' < 0, but to a 
much smaller extent at Xp > 0- The situation with unfavorable resonances (not containing 
the current quark) can be rather different. An interesting example is $ meson. 
Note first of all, that the neutrinoproduction of direct $ mesons, both unfavorable and espe- 
cially favorable, is expected to be strongly suppressed as compared to other favorable vector 
mesons, for example, u meson (with the same spin and isospin) which can be produced via 
the main subprocess ud followed by the recombination of the current u- quark with a 

u- antiquark created at the string breaking. The production of a favorable $ meson, which 
proceeds via the subprocess vu fi~s followed by the recombination with a s- quark, is 
triple-suppressed as compared to the u production due to: i) the suppressed magnitude of 
the Vus element of Cabibo-Kobayashi-Maskawa matrix [2]; ii) the strangeness suppression 
in the quark string fragmentation [3]; iii) the smallness of the fraction of the nucleon mo- 
mentum (integrated over the Bjorken xb variable) carried by the u- antiquark as compared 
to that for the d- quark (see e.g. [1]). The product of these suppression factors is about 
5 • 10^^. Even when taking into account the contribution from the competing recombination 
processes, e.g. uu n^, p^,ri,ri for the case of u and ss ^ r],!] for the case of 0, the 
expected ratio of the yields of favorable and u mesons remains very small, no more than 

lo-l 

This value is in strict contrast with experimentally estimated ratio {n^{xp > 0))/ {n^{xp > 
0)) = 0.18 ± 0.08 a.t xp > [1], indicating on the dominant role of indirect processes in the 
$ neutrinoproduction. One of them is, probably, related to the production of the charmed, 
strange Dj' meson via the subprocess on a strange sea quark, us — > fi~c, followed by the 
recombination of the charm quark with the sea s- antiquark from the nucleon remnant, 
(cs) — > D'^. Since here no creation of an extra (qq) pair is necessary, the meson can 
carry the overwhelming fraction z of the current c- quark energy and then transfer a signif- 
icant part of the latter to the $ meson in D'^ $X decays (which occur with a summary 
rate 16.1±1.6% [2]). An alternative subprocess, leading to the _D+ production (with ap- 
proximately the same probability as in the previous case), is ud — > fi~c followed by the 
recombination of the c- quark with a s- antiquark created at the string breaking. In this 



1 



case, however, Df carries on an average a lesser energy fraction z. 

This work is devoted to the experimental study of the neutrinoproduction of $ and Df 
mesons, with a particular aim to check the aforesaid mechanism of the indirect $ production. 
In Section 2, the experimental procedure is described. Section 3 presents the experimental 
data on the total and differential yields of $ meson. In Section 4, several decay modes of 
meson are analyzed and an estimation of its yields is inferred for the case of the leading 
meson carrying the overwhelming fraction {z > 0.85) of the current quark energy. The 
results are summarized in Section 5. 



2 Experimental procedure 

The experiment was performed with SKAT bubble chamber exposed to a wideband 
neutrino beam obtained with a 70 GeV primary protons from the Serpukhov accelerator. 
The chamber was filled with a propane-freon mixture containing 87 vol% propane (Csi^s) 
and 13 vol% freon (CF^Br) with the percentage of nuclei H:C:F:Br = 67.9:26.8:4.0:1.3 %. 
A 20 kG uniform magnetic field was provided within the operating chamber volume. 
Charged current (CC) interactions containing a negative muon with momentum >0.5 
GeV/c were selected. Other negatively charged particles were considered to be vr" mesons, 
except for the cases explained below. Protons with momentum below 0.6 GeV/c and a 
fraction of protons with momentum 0.6-0.85 GeV/c were identified by their stopping in 
the chamber. Non-identified positively charged particles were considered to be tt~^ mesons, 
except for the cases explained below. Events in which errors in measuring the momenta 
of all charged secondaries and photons were less than 60% and 100%, respectively, were 
selected. The mean relative error (Ap/p) in the momentum measurement for muons, pions 
and gammas was, respectively, 3%, 6.5% and 19%. Each event is given a weight which 
corrects for the fraction of events excluded due to improperly reconstruction. More details 
concerning the experimental procedure, in particular, the estimation of the neutrino energy 
El, and the reconstuction of 7r° — * 2'j and neutral strange particle decays can be found in 
our previous publications [SI d E]- 

The events with 3 < i?,^ < 30 GeV were accepted, provided that the reconstructed mass W 
of the hadronic system exceeds Wmin = 2 GeV. No restriction was imposed on the transfer 
momentum squared Q^. The number of accepted events was 4577 (5784 weighted events). 
The mean values of the kinematical variables were (E^) = 10.7 GeV, (W) = 3.0 GeV, (W"^) 
= 9.6 GeV2, (Q2) = 2.8 (GeV/c)^. 

About 8% of neutrino interactions occur on free hydrogen. This contribution was subtracted 
using the method described in [9] [TO] . 

When considering the production of resonances decaying into charged kaons, the K~ and 
hypothesis was applied, respectively, for negatively charged particles and non-identified 
positively charged particles with momenta p > Pcut = 0.55 GeV/c, introducing thereat 
proper corrections for the momentum of these particles. It has been checked that the choice 
of lower values of pcut does no practically influence the results presented in next sections. 
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3 The total and differential yields of (f) meson 



The K^K~ effective mass distributions for three different ranges of xp are plotted in Figure 
1. The main problem in the (j) signal separation is the large background from misidentified 
7r~^7T~ pairs having small effective masses. The background related to the correlated low- 
mass TT^vr" from rj — > 7r"'"7r~7r°, lj — > 7r"'"7r~7r° and p — > 7t~^71~ decays was subtracted from 
experimental distributions using the data on these resonances obtained in [1]. To describe 
the remaining (combinatorial) background, two methods were applied: 

i) The mixed event method, in which the shape of the background distribution was determined 
combining kaons from different events of the same topology and ranges of global kinematical 
variables {E,^, W). The normalization of the background distribution was kept as a free 
parameter when fitting the experimental mass distribution; 

ii) The fitted background method, in which the background was approximated by 

BG{m) = B ■ [m — 2mx)^ ■ exp (—7m) , (1) 

where rriK is the charged kaon mass, while B, (3 and 7 are free parameters (in some cases 7 
was fixed to zero). 

The 'signal' distribution was parametrized by a Gaussian form 

G^(m) = a^-exp[- ^^~7'^^' ] , (2) 

with the fixed pole mass = 1019 MeV and the experimental resolution cr^ = 4 MeV 
(estimated by simulations, resulting in a FWHM value F^^** ^ 10 MeV exceeding largely the 
(j) natural width F^ ^ 4.3 MeV). 

The mass distributions shown in Figure 1 were fitted as a sum of the 'signal' and background 
distributions. The fit results for both methods i) and ii) turned out to be practically the 
same. As it is seen, the data exhibit a clear signal aX xp > 0, with the corresponding 
yield {nq>{xF > 0)) = (0.92 ± 0.34) ■ 10~^ (corrected for unobserved decay modes). No 
production is observed in the backward hemisphere {xp < 0), where the data are consistent 
with the background distribution. As a result, the signal in the full xp- range turns out 
to be less expressed as compared to that at xp > 0, leading to a less accurate estimation of 
the total yield: = (1.19 ± 0.61) ■ lO'^. 

We have also obtained the Wmin- dependence of {nq,{xp > 0)) and found it to increase 
continuously from {n^{xp > 0)) = (0.8 ± 0.3) ■ 10~^ at Wmin = 1-8 GeV [Ij to {n^{xp > 
0)) = (1.23±0.53)-10-2 and (1.44±0.69) ■ 10"^ at W^in = 2.6 GeV and 2.9 GeV, respectively. 
The differential spectrum dn^/dxp, compared to that for p° meson [11] (at > 2 GeV), 
is presented in Figure 2, from which two observations emerge. Firstly, unlike $ meson, 
the yield of in the backward hemisphere is not much smaller than that in the forward 
hemisphere. This fact can be explained [llj by a contribution from secondary intranuclear 
interactions of pions, vrA^ — > p°X, resulting in a p° multiplicity gain, mainly at xp < 0. 
Nothing similar happens with 0, due to a much smaller probability of the production 
in secondary interactions (see Section 5 for details). Secondly, the yield of p° is rather 
prominent at large xp > 0.6 — a characteristic feature for the direct leptoproduction of a 
favorable hadron. On the contrary, the yield of is more concentrated around moderate 
positive values of Xi;', indicating on its origin from the decay of a favorable leading resonance. 
A probably candidate for the latter is, as it was mentioned in Introduction, the charmed, 
strange meson the production of which is considered in the next section. 
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4 The yield of meson 



As it was mentioned above, the xp- dependence of the meson yield indicates on its origin 
from the decay of leading mesons which can be produced, for example, in exclusive 
reactions 

pN ^ li- + N* + Dj {or , (3) 

and 

z/iV ^ /i- + r + D+ {or D*+) , (4) 

with the vector D*"^ state decaying into D'^j (94.2%) or D^n^ (5.8%), a nucleon remnant 
N* (turning into a nucleon or a low-state resonance), and a hyperon Y . 
In reactions (4) and (5), most of D'^ mesons carry the overwhelming fraction z of the current 
c- quark energy. The mean yield of the leading vector meson (with z > 0.75) was mea- 
sured earlier in uNe- interactions, resulting in {nj-,*-{z > 0.75)) = 0.051 ± 0.016 One 
can assume, that the mean yield of the pseudoscalar D~ (direct production) composes 1/3 
of that for D*~, i.e. the summary yield of is expected to be about 0.068 ± 0.021, leading 
to the mean yield of the decay (j) mesons {n'^'^) = 0.014 ± 0.006 per CC vNe- interaction. 
Below we consider several hadronic decay modes of D"^ with branching ratios exceeding a 
few percent (see Table). A rather severe cut was applied on z {> 0.85) in order to reduce 
the contribution from the background processes. The events with W < 2.9 GeV, i.e. below 
the threshold of the production, were excluded (it should be, however, noted, that the 
experimental resolution of W is about 10% in this experiment, and hence the events with 
the estimated value of W in between 2.6 < W < 2.9 GeV can also have a small contribution 
to the D'^ production). The experimental mass resolutions (estimated by simulations) for 
considered channels are quoted in Table. 

Table. The considered D'^ decay modes, experimental mass resolutions, estimated yields 
at z > 0.85 and decay fractions. 



Decay mode 


Mass resolution (MeV) 


Estimated yield 


Decay fraction 






(in 10-2) 


(in %) [2J 




22 


0.56±0.26 


5.50±0.28 


b) K+K-TT+TT° 


45 


0.50±0.25 


5.6±0.5 


c) 37r+27r-7r° 


42 


0.17±0.12 


4.9±3.2 


d) p^rj', rj' p°7 


52 


0.25±0.20 


3.6±0.6 


e) p^rj^rj -k^-k^-k^ 


120 


0.09±0.09 


3.0±0.5 


f) p+r/,r/ 77 


200 


0.27±0.20 


5.1±0.8 


sum 




1.84±0.48 


27.7±3.4 



The main channels with a K^K'^ pair in the final states are: a) K^K't^^ and b) 
K^K'n^TT^. The contamination from tt+tt" pairs is expected pt2] to be more prominent at 
large (in the absolute value) negative values of cos??|^^, "&*kk being the angle between the 
K^K~ direction in the D'^ candidate rest frame and the direction of the Lorentz boost from 
the lab system to the rest frame. We excluded the combination with cos'd*^^^ < —0.9 



4 



where a significant enhancement of the number of combinations was observed. For channel 
b) containing a 7r°, a cut 115 < m^^ < 155 MeV was apphed for the 77 effective mass (the 
same cut was also applied for channels c)— f) containing one or two neutral pions in the final 
state). The effective mass distributions for channels a) and b) are plotted in Figure 3. For 
both channels, an excess of events around the D'^ mass is observed above the combinatorial 
background which can be satisfactorily described using the mixed event method (the dashed 
curves). The experimental distributions were fitted by a sum of a Gaussian (with the width 
quoted in Table) and the mixed background distribution. The resulting yields, corrected for 
the applied cuts and the contamination from the background 77 combinations (see [13] for 
details), are quoted in Table. 

A rough estimation for the yield of the channel c) was extracted from the excess of events 
above the combinatorial background in the mass region (Figure 3 c). The combinatorial 
background was obtained with the mixed event method and normalized to the experimental 
distribution above 2.15 GeV/c^ (dashed curve). The mean yield, obtained after introducing 
corrections related to the 7r° reconstruction, is quoted in Table. 

To select events-candidates to channel d) p^ri'{rj' — > p°7), the following cuts were applied: 
the effective mass of tt^tt^ from the p° decay is enclosed in the range 0.6—0.95 GeV/c^; the 
effective mass of tt+tt^ from the decay is enclosed in a slightly wider (as compared to p°) 
range of 0.57—0.97 GeV/c^, in view of worse mass resolution (crp+ = 47 MeV) as compared 
to (TpO = 20 MeV (see ^ for details); the effective mass of p°7 from the rj' decay is enclosed 
in the range of m^/ ± 2(j{ri' p^l), with the mass resolution a{r]' —>■ p^'-f) = 30 MeV. The 
effective mass distribution for channel d) plotted in Figure 3 d is fitted as a sum of a Gaus- 
sian (with the width quoted in Table) and the mixed background distribution. The resulting 
yield, corrected for the applied cuts, is quoted in Table. 

To select events-candidates to channels e) p^rj{r] — > 7r"'"7r~7r°) and f ) p^f]{ri 77) the follow- 
ing cuts were applied: the effective masses of Tr+vr^vr" and 77 from the 77 decay are enclosed 
in the ranges of ± 2a{ri tt+tt^tt^) and ± 2a{ri 77), with the mass resolutions, 
respectively, a{r] — ^ tt+tt^tt^) = 28 MeV [Ij and a{ri — > 77) = 62 MeV. After these cuts, only 
one event-candidate to channel e) and two event-candidates to channel f ) survived. All three 
events had effective masses compatible with the D'^ mass within ±1.5 mass resolution for 
the corresponding channel (see Table). The corresponding yields, corrected for the applied 
cuts are quoted in Table. 

Finally, we looked for decay channels with a in the final state: K~^Kg, KgK~7r~^7r~ and 
K+K'^TT+TT-, with the summary branching ratio (4.09±0.22) % [2]. No combination with 
z > 0.85 was found in the mass region, probably due to smallness of their branching 
ratios and the restricted statistics of mesons in our experiment. 

The summary experimental yield of channels a)— f) (with the summary branching ratio 
(27.7±3.4) %) is equal to (1.84 ± 0.48) ■ lO'^, resulting to the total yield of leading D+ 
mesons equal to {nj^+{z > 0.85); W > 2.9 GeV ) = (6.64 ± 1.91) • 10"^. This value can be 
compared with the mean yield of the vector D*~ in i?Ne- interactions measured at slightly 
different kinematics: (n^.-(z > 0.85); W >3 GeV ) = (5.1 ± 1.6) ■ lO'^ [12]. 

5 Discussion and Summary 

The inclusive production of (j) mesons in neutrinonuclear interactions is studied for the first 
time. It is found that the production occurs predominantly in the forward hemisphere of 
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the hadronic c.m.s., with the mean yield varying from {n^{xF > 0)) = (0.92 ± 0.34) ■ 10 ^ 
at W > 2 GeV up to (1.23 ± 0.53) ■ 10"^ at W > 2.6 GeV and (1.44 ± 0.69) ■ lO^^ at 
W > 2.9 GeV. The measured yields are much larger than expected for the case of the direct 
neutrinoproduction, hence indicating on a dominant role of indirect mechanisms. 
A possible candidate to the latter are the secondary intranuclear interactions like 

M + N ^(j) + X, (5) 

where M denotes pions, kaons, non-strange {p,uj) and strange (i^*(892)) vector mesons, 
as well as r] and rj' mesons. To estimate the multiplicity gain {n{M — >■ 0)) of meson 
in reactions (5), we used a simple model described and applied in [TOl [TTl [H]. The main 
quantities determining (n(M — > 0)) are the mean multiplicities {uMixp > 0)) of intermediate 
mesons produced in uN- interactions at Xi;' > and the cross sections a{M 0) of reactions 
(5). 

The values of {nT^{xp > 0)) were taken from [10], subtracting the contribution from the 
decays of mesonic resonances [Ij. The values of a{7r —>■ 0) were taken from [15j and averaged 
over the pion momentum range from ~ 2 GeV/c up to ~ 12 GeV/c (above which the 
pion yield is negligible in this experiment). The mean value of a{7i —>■ 0) averaged over the 
pion and target nucleon species is in the range of a^n ^ 0) = 0.035 ± 0.015 mb, resulting 
in {n{n 0)) = (0.05 ± 0.03) ■ 10~^. The cross sections a{p 0) and a{u! 0) are 
expected to be slightly larger as compared to a{7r — >■ 0) due to larger mean momenta of p 
and LV. Nevertheless, due to the smallness of their mean multiplicities [1], the contribution 
to the production is approximately by one order of the magnitude smaller as compared to 

{n{7T 0)). 

The averaged cross section for kaon-induced reactions, estimated from the data compiled in 
|16j . is in the range of a{K — >■ 0) = 0.12 ± 0.04 mb. The mean multiplicity of charged kaons 
is assumed to be equal to that of neutral ones taken from [14J. The mean multiplicities of 
vector mesons i^'*(890) were taken from [1], while the cross section a{K* —>■ 0) was assumed 
to be equal to a{K — 0). As a result, the summary contribution {n{K — > 0)) + {n{K* — > 0)) 
was estimated to be (0.01 ± 0.01) ■ 10"^ 

The mean multiplicities of t] and 1]' were taken from [Ij. In view of the ss content of rj and rj' 
(with probabilities 2/3 and 1/3, respectively), the production cross sections were assumed 
to be related to that of p° production by pions: a{ri — > 0) ;^ 2/3(T(7r ^ p^) ^ 2 mb and 
a-{ri' — > 0) ^ l/3(T(7r — > p°) ^ 1 mb (at p-,, around a few GeV/c). Under this assumption, 
one obtains {n{ri 0)) = (0.102 ± 0.054) ■ lO'^ and (n(r/' ^ 0)) = (0.016 ± 0.010) ■ lO'^, 
where the uncertainty in cross sections is not included in the quoted errors. The summary 
contribution from all aforesaid processes does not exceed 0.2 ■ 10~^. This value, even when 
assuming that all produced mesons acquire xp > (which is not the case), composes 
only a minor fraction of the measured yield {n^i^xp > 0)) = (0.92 ± 0.34) ■ 10~^ a.t W > 2 
GeV. An appreciable contribution to the latter can be provided if only the unknown cross 
sections a{ri — > 0) and a{ri' 0) were a few times larger than supposed above. It might be 
worthwhile to note, that an information about the role of secondary intranuclear interactions 
can be, in principle, inferred from the data (not yet existing) on the A- dependence of the 
yield in neutrinonuclear reactions. 

Another source of indirect production at xp > are the decays of leading mesons 
carrying the overwhelming fraction z of the hadronic energy. To look the neutrinoproduction 
of D'^ mesons, six different final states of decays are considered requiring z > 0.85. For 
all of them, an excess of events is observed at the mass region, allowing to infer, for the 
first time in neutrino-induced reactions, an estimation for its yield: {D'^{z > 0.85); W > 2.9 
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GeV ) = (6.64 ±1.91) ■ IQ-^. 

The events-candidates of the production were collected to a subsample used furthermore 
for simulation of the ^ (p + and Df ^ + vr"'" decays composing the overwhelming 
fraction (85%) of the D'^ (pX decays The simulated differential spectrum dn'^p'^/dxp 
turned out to be strongly shifted towards the forward hemisphere in the hadronic c.m.s., 
providing 85.4% of mesons to have xp > ^ and resulting in the expected mean yield 
of decay mesons equal to {nf^xp > 0);W > 2.9 GeV ) = (0.91 ± 0.26) ■ 10"^. This 
composes more than (or about) the half of the measured value {n^^xp > 0)]W > 2.9 
GeV ) = (1.44 ±0.69) ■ lO'^. 

In order to make a comparison with the measured differential spectrum at > 2 GeV 
(plotted in Figure 2), the yield of D'^ was properly renormalized to the number of events 
with W > 2 GeV, assuming that the D'^ production at W < 2.9 GeV can be neglected. This 
lead to {n^+{z > 0.85); W >2 GeV ) = (2.74 ± 0.79) ■ 10"^ resulting in the expected yields 
of the decay mesons (n^"'=(all xp);W > 2 GeV ) = (0.44 ± 0.13) ■ 10"^ and {nf%xp > 
0);W > 2 GeV ) = (0.38 ± 0.11) ■ 10"^. These compose less than (or about) the half of 
experimentally estimated values, respectively, (n(^(all xp);W > 2 GeV) = (1.19±0.61)-10^^ 
and {n^{xp > 0);W > 2 GeV ) = (0.92 ± 0.34) ■ lO^^. As a result, the magnitude of the 
predicted spectrum dn'^'^/dxp approximately twice underestimates that for the measured 
one (Figure 2). The shape of the predicted spectrum is, however, compatible with the 
experimental one. 
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Figure 1: The K^K~ effective mass distributions aX W > 2 GeV. The curves arc the fit 
results for two cases of the background description: using the analytical form (1) (the left 
panel) and the mixed event method (the right panel). 
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Figure 2: The xp spectra for cf) and p° at > 2 GeV. The curve is the expected distribution 
for (f) mesons from the Df decays. 
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Figure 3: The effective mass distribution for channels a)— d) at W > 2.9 GeV. The curves 
are the fit results (see text). 
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